Confinement- Induced Resonances in Low-Dimensional Quantum Systems 



Elmar Haller, 1 Manfred J. Mark, 1 Russell Hart, 1 Johann G. Danzl, 1 Lukas 

Reichsollner, 1 Vladimir Melezhik, 2 Peter Schmelcher, 3 and Hanns-Christoph Nagerl 1 

1 Institut fiir Experimentalphysik and Zentrum fiir Quantenphysik, 
Universitdt Innsbruck, Technikerstrafie 25, 6020 Innsbruck, Austria 
2 Bogoliubov Laboratory of Theoretical Physics, Joint Institute for Nuclear Research, Dubna, 141980 Dubna, Russia 
Zentrum fiir Optische Quantentechnologien, Universitiit Hamburg, Lumper Chaussee 149, 22761 Hamburg, Germany 

(Dated: April 6, 2010) 

We report on the observation of confinement-induced resonances in strongly interacting quantum-gas systems 
with tunable interactions for one- and two-dimensional geometry. Atom-atom scattering is substantially modi- 
fied when the s-wave scattering length approaches the length scale associated with the tight transversal confine- 
ment, leading to characteristic loss and heating signatures. Upon introducing an anisotropy for the transversal 
confinement we observe a splitting of the confinement-induced resonance. With increasing anisotropy additional 
resonances appear. In the limit of a two-dimensional system we find that one resonance persists. 
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Low-dimensional systems have recently become exper- 
imentally accessible in the context of ultracold quantum 
gases. For a two-dimensional (2D) geometry, the Berezinskii- 
Kosterlitz-Thouless (BKT) transition has been observed (TJ, 
and in one dimension the strongly-correlated Tonks-Girardeau 
(TG) [2-6] and super- Tonks-Girardeau (sTG) gases [6| have 
been realized. In these experiments steep optical potentials 
freeze out particle motion along one or two directions and re- 
strict the dynamics to a plane or to a line. Such quasi-2D or 
quasi- ID systems can be realized with ultracold gases when 
the kinetic and the interaction energy of the particles are insuf- 
ficient to transfer the particles to transversally excited energy 
levels. Whereas the confinement removes motional degrees 
of freedom, it also provides an additional structure of discrete 
energy levels that can be used to modify scattering along the 
unconfined direction and by this to effectively control the in- 
teraction properties of the low-dimensional system ||7}j9). In 
this Letter, we investigate the few-body scattering processes 
that give rise to the capability to tune interactions and hence 
to drastically alter the properties of low-dimensional many- 
body quantum systems [6|. 

In three-dimensional (3D) geometry magnetically-induced 
Feshbach resonances (FBRs) [10] allow tuning of the inter- 
particle interaction strength. A FBR occurs when the scat- 
tering state of two atoms is allowed to couple to a bound 
molecular state. Typically, scattering state and bound state are 
brought into degeneracy by means of the magnetically tun- 
able Zeeman interactions. For particles in ID and 2D geom- 
etry a novel type of scattering resonance occurs. Coupling 
between the incident channel of two incoming particles and 
a transversally excited molecular bound state generates a so- 
called confinement-induced resonance (CIR) 1714^ fTTI4T3l . A 
CIR occurs when the 3D scattering length approaches the 
length scale that characterizes the transversal confinement, i.e. 
the harmonic oscillator length «x = uhj (mco±) for a parti- 
cle with mass m and transversal trapping frequency co±. This 

causes the ID coupling parameter gw = 2 ^" 2 3D \-cl Ja to 



diverge at a± — Cam, where C = 1.0326 is a constant (TIED. 
The CIR allows tuning of interactions from strongly repulsive 
to strongly attractive and thus represents a crucial ingredient 
for the control of interactions in a low-dimensional system. 
Modification of scattering properties due to confinement has 
been measured near a FBR for fermions 1 14 1, and, recently, a 
CIR has been observed for a strongly-interacting ID quantum 
gas of bosonic Cs atoms and was used to drive the crossover 
from a TG gas with strongly repulsive interactions to an sTG 
gas with strongly attractive interactions |6j- Here, for an ul- 
tracold quantum gas of Cs atoms with tunable interactions, we 
study the properties of CIRs by measuring particle loss and 
heating rate and, in particular, confirm the resonance condi- 
tion a± = Ca3D for symmetric ID confinement. For the case 
of transversally anisotropic confinement we find that the CIR 
splits and, to our surprise, persists for positive 030 even when 
the anisotropy reaches the limit of a 2D system. 

Figure [TJa) reviews the basic mechanism that causes a CIR 
for zero collisional energy in ID (9). It is assumed that in 
3D the scattering potential supports a single universal bound 
state for strong repulsive interactions (dotted line) iflOl . The 
point where the incoming channel of two colliding atoms and 
the universal dimer state are degenerate marks the position 
of a 3D FBR (triangle). In ID, strong transversal confine- 
ment shifts the zero-energy of the incoming channel (middle 
dashed line) and introduces a transversally excited state (up- 
per dashed line). As a result of the strong confinement, which 
modifies the long-range part of the molecular potential, the 
universal dimer state with binding energy £3 (lower solid line) 
exists also for attractive interactions [15] whereas the original 
3D FBR has disappeared. Instead, there is a CIR (star) when 
the incoming scattering channel becomes degenerate with the 
transversally excited molecular bound state (upper solid line). 
It is assumed that the binding energy of this state is also E%, 
shifted by 2ha>±\7\. In more detail, as depicted in Fig.[TJb), 
we assume that the energy levels of non-interacting atoms, as a 
result of cylindrically symmetric transversal confinement, can 
be approximated by those of a 2D harmonic oscillator with 
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FIG. 1: (color online) (a) Illustration of the mechanism responsible 
for a CIR, see Ref . 1 9 1 and text for details. The energy levels near 
a scattering resonance are plotted as a function of 1/«3d- The CIR 
occurs for C«3d = a±_ when scattering atoms are allowed to cou- 
ple to transversally excited bound states, (b) indicates the shift and 
splitting for anisotropic confinement characterized by Aco. (c) Ex- 
perimental configuration. Two laser beams create an optical lattice 
that confines the atoms to an array of approximately 3000 indepen- 
dent, horizontally-oriented elongated ID tubes, (d) Tuning of 030 is 
achieved by means of a FBR with a pole at B = 47.78(1) G fl8l . 



Em ,m = ^ (°± ( n 1 + n 2 + 1 ) and quantum numbers n \ and n.2 be- 
longing to the two Cartesian directions. Scattering atoms ifTBI 
in the transversal ground state (0,0) can couple to the excited 
states («i,«2) if the parity of the total wave function is pre- 
served [12|. The energetically lowest allowed excited states 
are threefold degenerate with an energy E = 3hco± and with 
quantum numbers (1,1), (2,0) and (0,2). For the transver- 
sally symmetric confinement, they contribute towards a single 
CIR (9). However, the contribution of the state (1,1) is neg- 
ligible due to the zero contact probability of the atoms and 
the short-range character of the interatomic interaction. Un- 
equal transversal trapping frequencies C0i and C02 = C0\ + Aco 
lift this degeneracy and shift the energy levels according to 
E ni .„ 2 — ftG>i ( n i + n 2 + 1) +hAco(ri2 + 1 /2). One thus expects 
a splitting of the CIR. 

We start from a tunable Bose-Einstein condensate (BEC) 
of 1 .0 to 1 .4 x 10 s Cs atoms in the energetically lowest hyper- 
fine sublevel |[T7l confined in a crossed-beam optical dipole 
trap and levitated against gravity by a magnetic field gradi- 
ent of |VZ?| ss 31.1 G/cm. Tunability of is given by a 
FBR as shown in Fig. [TJd) with its pole at Bq = 47.78(1) 
G and a width of 164 mG ifTTl [TSI . The resonance resides 
on top of a slowly varying background that allows tuning of 
fl3D from to values of about 1000 an, where an is Bohr's 
radius. Using the FBR, we can further tune to values up 
to 030 ~ 6000 an given our magnetic field control with an 
uncertainty of AB « 10 mG. We convert B into ajo using 
the FBR parameters from Ref. [ 1 8 1 . The BEC is produced 
at C3D ~ 290 an. We load the atoms within 300 ms into an 
optical lattice, which is formed by two retro-reflected laser 
beams at a wavelength of A = 1064. 49( 1 ) nm and with a beam 



FIG. 2: (color online) Particle loss and heating rates in the vicin- 
ity of a CIR. (a) The number N of remaining atoms after T = 200 
ms shows a distinct drop ("edge") when B is scanned across the 
CIR. A clear shift of the position of the edge to lower values for 
B can be observed when the transversal confinement is stiffened, 
<a ± = 2% x (0.84,0.95, 1.05) x 14.2(2) kHz (circles, squares, trian- 
gles), (b) Position of the edge (circles) as determined from the in- 
tersection point of a second-order polynomial fit to the minimum for 
N and the initial horizontal baseline as shown in (a), converted into 
values for a^r,. For comparison, the position of the minimum (trian- 
gles) is also shown. The solid line is given by Ca^u = a± with the 
predicted value C = 1 .0326. (c) Heating rates near the CIR (circles). 
For comparison, N is also shown (triangles). For this measurement, 
tt>x — 2tt x 12.0(2) kHz. All error bars reflect lcr statistical uncer- 
tainty. 



waist of approximately 350 Lira, one propagating vertically 
and one propagating horizontally as illustrated in Fig. [TJc). 
These lattice beams confine the atoms to an array of approxi- 
mately 3000 horizontally oriented, elongated ID tubes with a 
maximum occupation of 60 atoms at a linear peak density of 
approximately ps 2/jj.m. Weak longitudinal confinement 
results from the Gaussian-shaped intensity distribution of the 
beams. We raise the lattice to a depth of typically V = 30 Er, 
where En = h 2 / (2mA 2 ) is the photon recoil energy. At this 
depth, the resulting transversal and longitudinal trap frequen- 
cies are co± = 2% x 14.5 kHz and 0)|| = 2n x 16 Hz and we 
then have a± sa 1370 an. After loading we slowly ramp down 
|VB| in 50 ms and adiabatically increase a3D to 915 an in 100 
ms to create a TG gas with well-defined starting conditions 
near the CIR |6|. To detect the CIR as a function of B, mani- 
fested by a loss resonance, we quickly set B in less than 200 lis 
to the desired value, wait for a hold time of typically z = 200 
ms, and then measure the number N of remaining atoms by 
absorption imaging. For this, we re-levitate the atoms, ramp 
down the lattice beams adiabatically with respect to the lattice 
band structure, and allow for 50 ms of levitated expansion and 
2 ms time-of-flight. Note that T is chosen to be much longer 
than the lifetime of the sTG phase J6) . 

We observe the CIR in the form of an atomic loss signature 
as shown in Fig. [2] We attribute the loss near the resonance 
to inelastic three -body |fl9l or higher-order collisions ll20l . 
which lead to molecule formation and convert binding energy 
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into kinetic energy, causing trap loss and heating, similar to 
the processes observed near a FBR [ 10 1 . In our case, inelastic 
two-body processes can be ruled out for energetic reasons and 
single-particle loss occurs on the timescale of tens of seconds. 
In Fig. |2ja) the CIR can be identified as a distinct "edge" for 
the atom number N. Initially, in the TG regime of strong re- 
pulsive interactions, here for B < 47.35 G, losses are greatly 
suppressed, but increase rapidly on the attractive side of the 
CIR. N drops to a minimum when B is increased and then re- 
covers somewhat. A clear shift of the loss signature to lower 
values for B and hence lower values for can be discerned 
when the confinement is stiffened. When we identify the po- 
sition of the edge with the position of the CIR, we find good 
agreement with the analytical result Co3d = a± as shown in 
Fig.[2jb). As we have no theoretical description of the detailed 
shape of the loss resonance, we also plot, for comparison, the 
position of the minimum, which is shifted accordingly. 

In Fig. |2jc) we juxtapose the loss and the heating rate that 
we measure in the vicinity of the CIR. For this, we measure 
the increase of the release energy within the first 100 ms. Af- 
ter holding the atoms for time T at a given value of B, we 
decrease back to 250 an in 20 ms, switch off the lattice 
potential and determine the release energy in the direction of 
the tubes from the momentum distribution in free space ex- 
pansion. We observe an increase for the heating rate when 
the CIR is crossed. From a low value of 10 nKVs in the TG 
regime it rises to a maximum of approximately 150 nK/s and 
then drops to settle at some intermediate value. The position 
of the maximum agrees well with the maximum for atom loss. 
We check that the system's increase in energy is sufficiently 
small so that its ID character is not lost. The release energy, 
even at maximal heating, remains below kg x 30 nK, which is 
far below the energy spacing of the harmonic oscillator levels, 
ha>± ~ijx 600 nK. 

We now examine ID systems with transversally anisotropic 
confinement. Starting from a lattice depth of V = 25 Er 
along both transversal directions, yielding co± = ©i = ©2 = 
2n x 13.2(2) kHz, we increase the horizontal confinement to 
frequencies up to ©2 = 27T x 16.5(2) kHz, corresponding to a 
lattice depth of 39 Er, while keeping the depth of the verti- 
cal confinement constant. Fig. |3ja) shows a distinct splitting 
of the original CIR into two loss resonances, CIRi and CIR2. 
The splitting increases as the anisotropy is raised. In Fig. [3jb) 
we plot the 3D scattering length values 030,1 an d a 3D,2 that we 
associate with the positions of CIRi and CIR2 as a function of 
the frequency ratio (£h/(i>\. For this, as it becomes difficult 
to assign an edge to both of them, we simply determine the 
associated atom number minima and subtract a constant off- 
set of 88(7) fl() as determined from the measurement shown in 
Fig.|2jb). One of the resonances, CIR2, exhibits a pronounced 
shift to smaller values for 030 as the horizontal confinement is 
stiffened. The second resonance, CIRi, shows a slight shift 
towards higher values for ajo- We now use the lifting of 
the degeneracy for the energy levels as indicated in Fig. |TJb) 
to model the observed splitting of the CIR. We assume that 
the implicit equation £(1/2,— £B/(2ft©x) + l/2) = — flx/ a 3D 



x 10 4 



10 



-Q 

E 



10 



-Q 

E 



(a) 



CIR "V* 



10 z 

8 a 

6 I 
c 

4 E 
o 

2 ™ 



47.1 



47.2 



47.3 



47.4 



47.5 



47.6 



47.7 47.8 



magnetic field B (G) 




1.00 1.05 1.10 1.15 1.20 1.25 0.0 
frequency ratio (a/to, 



0.1 0.2 

Aco/to, 



0.3 



FIG. 3: (color online) Splitting of a CIR for a ID system with 
transversally anisotropic confinement, (a) As the horizontal confine- 
ment is stiffened, a^/cQi = 1.00, 1.10, 1.18 (circles, diamonds, trian- 
gles) for (Q{ = 1% x 13.2(2) kHz, the CIR splits into CIRi and CIR 2 . 
(b) Position of CIRi («3d,i, circles) and CIR2 (fl3D,2, squares) as a 
function of the frequency ratio (Oi/(0\. (c) Binding energy difference 
A£b as determined from the implicit equation (see text) in compar- 
ison to the expectation from the simple harmonic oscillator model 
(solid line). 



for the binding energy E% [9 | remains approximately valid for 
sufficiently small A©, taking © ± = d>\- Here, £ is the Hur- 
witz zeta function. We translate the scattering length values 
a3D,i and «3d,2 into binding energies and calculate the energy 
difference A£b = £b(«3D.i) — E^{a^\)^), shown in Fig. |3jc). 
While this model does not explain the upward deviation seen 
for CIRi, the difference A£b is in reasonable agreement with 
the expected energy shift caused by the shifts of the excited 
harmonic oscillator states (£0.2 —£2.0) = 2hAco (solid line in 
Fig. [3jc)). We thus attribute CIR2 to the stiffened confine- 
ment along the horizontal direction and hence to state (0,2), 
while CIRi corresponds to the unchanged confinement along 
the vertical direction and hence to state (2,0). 

We observe the appearance of additional structure in 
the measured loss curves when we increase the transversal 
anisotropy by weakening the confinement along one axis, here 
along the vertical direction. Fig.Qa) shows the atom number 
after 1 = 300 ms for trapping frequency ratios ©1 /©2 from 
0.67 to 0.45. Multiple loss resonances appear close to the po- 
sition of CIRi. The number of resonances increases and the 
positions shift continuously as the confinement is weakened. 
We speculate that those resonances are a result of a coupling 
to additional excited states, resulting in a multi-channel scat- 
tering situation. Also the weakening of the confinement could 
induce sufficient anharmonicity to allow for violation of the 
parity rule BTl . 
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FIG. 4: (color online) (a) Appearance of additional structure in the 
vicinity of CIRs for strongly anisotropic transversal confinement. 
The trap frequencies are a>2 = 2k x 16.6(2) kHz and a>\/o>i = 
0.67,0.60,0.53,0.49,0.45 from top to bottom, (b) Transition from 
ID to 2D confinement. As the horizontal lattice is ramped down, 
CIR2 shifts and persists, while CIRj disappears (C0\ = 2n x 13.0(2) 
kHz and ©2/0)] = 0.58, 0.42, 0.00 for squares, circles, and triangles), 
(c) Scaling of the CIR's position in 2D, in analogy to the ID case 
shown in Fig.[2jb). The position of the CIR as determined from the 
edge (circles) and, alternatively, from the minimum in atom number 
(triangles) shifts to lower values for a^u as the confinement is stiff- 
ened and flx,2D is reduced. The solid line is a linear fit according to 
C2D03D = «±,2D with C 2D = 1.19(3). 

Surprisingly, we find that one of the CIRs persists in the 
limit of a 2D system. Previous theoretical studies on 2D sys- 
tems have predicted the appearance of a CIR for negative am, 
but not for positive 030 [ 13 , 22 1. In the experiment, we reduce 
the horizontal confinement while keeping the vertical confine- 
ment constant to probe the transition from the array of tubes to 
a stack of pancake-shaped, horizontally-oriented 2D systems. 
Trapping in the horizontal direction is still assured, now by the 
Gaussian profile of the vertically propagating laser beam, for 
which ©2 = 2tz x 11 Hz. Fig. |4}b) shows that the CIR asso- 
ciated with the tight confinement shifts to lower values for B 
and hence for 030 as the horizontal confinement is weakened. 
In the limit of 2D confinement, one of the CIRs, and in fact 
all the additional structure observed above, have disappeared, 
but one resonance persists. To check that the observed reso- 
nance is indeed the result of the 2D confinement, we vary the 
confinement along the tight vertical direction. Fig.|4|c) plots 
the positions of edge and minimum of the loss signature as a 
function of ax,2Di the confinement length associated with this 
direction. When we again associate the edge with the pole of 
the resonance, we obtain C2d«3d = fl ±.2D with C2D = 1-19(3), 
where C2D is a scaling factor similar to C for the ID case. Fur- 
ther scattering experiments are needed to elucidate the energy 
dependence of this 2D scattering resonance. 

In summary, we have investigated the properties of CIRs, 
which appear in low-dimensional quantum systems as a re- 



sult of tight confinement and which replace "conventional" 
3D Feshbach resonances to tune the effective atomic inter- 
action strength. We observed a splitting of the CIR for 
anisotropic transversal confinement, the appearance of mul- 
tiple resonances for strongly anisotropic confinement, and the 
survival of one resonance for positive am in the limit of 2D 
confinement. We expect that CIRs will not only be used in 
ID geometry to tune the effective interaction strength as re- 
cently demonstrated [6], but also in 2D geometry and mixed 
dimensions [23 1 for the study of strongly-interacting quantum 
systems. 
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